Observation of Multiple Thresholds in the Many- Atom Cavity QED Microlaser 
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We report the observation of multiple laser thresholds in the many-atom cavity QED microlaser. 
Traveling-wave coupling and a supersonic atom beam are used to create a well-defined atom-cavity 
interaction. Multiple thresholds are observed as jumps in photon number due to oscillatory gain. 
Although the number of intra-cavity atoms is large, up to N ~ f 3 , the dynamics of the microlaser 
agree with a single atom theory. This agreement is supported by quantum trajectory simulations 
of a many-atom microlaser and a semiclassical microlaser theory. We discuss the relation of the 
microlaser with the micromaser and conventional lasers. 



The most fundamental model of light-matter interac- 
tion at the atomic level consists of a two-level atom cou- 
pled to a single mode of the electromagnetic field, e. g. 
in an optical resonator. For an atom-cavity coupling 
strength greater than the decay rates of atom and cavity 
(.9 ^ r c ,r a ), an excited atom exchanges energy coher- 
ently with the cavity (Rabi oscillation). Q] 

Atom-resonator interaction is also the domain of laser 
physics; therefore it may be somewhat surprising that 
most descriptions of laser operation use an incoherent 
model involving population densities and Einstein A and 
B coefficients. Such a model applies due to gain medium 
broadening, laser field nonuniformity, and other statisti- 
cal effects 2] . 

The microlaser is to our knowledge the first laser in 
which coherent Rabi oscillation is explicitly reflected in 
the dynamics of the laser. A controlled atom-cavity inter- 
action and absence of strong statistical averaging leads to 
behavior foreign to conventional lasers. In this paper we 
report the most dramatic of these effects, the existence 
of multiple laser thresholds. 

The microlaser is the optical analogue of the micro- 
maser Q, in which a similar bistable behavior has been 
observed for a single atom in the cavity Q . The most sig- 
nificant differences between the two experiments are: (i) 
a large number of atoms is present, (ii) optical frequen- 
cies allow direct detection of generated light, and (iii) 
the number of thermal photons at optical frequencies is 
negligible. 

An earlier version of the microlaser experiment H was 
shown to exhibit laser action with an intracavity atom 
number N on the order of f . In the current setup, N 3> f ; 
remarkably, the many-atom system exhibits very similar 
behavior to that predicted from a single atom theory. 
This is supported by quantum trajectory simulations and 
may be explained by a semiclassical theory (0, 0). 

The experiment is illustrated in Fig. ^ A supersonic 
beam of 138 Ba atoms passes through the TEMoo mode 
of a high-finesse optical cavity (symmetric near-planar 
cavity, mirror separation L as 0.94 mm, mirror radius of 
curvature ro = 10 cm, finesse F as 9.0 x 10 5 .) The mirror 
separation and finesse were determined by measurement 
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FIG. 1: Schematic of cavity QED microlaser. M: cavity mir- 
rors, C: cavity mode, AB: atomic beam from oven, A: col- 
limating aperture, P: pump field (into page), 6: cavity tilt 
angle, LP: locking probe beam. Dashed line is normal to cav- 
ity axis. 



of transverse mode spacing and cavity ringdown decay 
time. The background pressure in the vacuum chamber 
is less than 10 -6 torr. Prior to entering the cavity mode, 
each atom is excited by a cw pump laser (Coherent 899- 
21 ThSapphire ring laser) locked to the 1 So <-> 3 Pi transi- 
tion (A = 791.1 nm, linewidth T a as 50 kHz) as described 
in . The pump beam is focused by a cylindrical lens to 
a roughly 50/im x 500/im elliptical Gaussian beam cen- 
tered approximately 150 /im from the cavity axis. This 
distance between pump and cavity mode assures that the 
pump beam does not overlap the cavity, while minimiz- 
ing the effect of atomic decay between the pump and 
cavity. The total decay rate from the 3 Pi excited state 
corresponds to a decay length as 1 .3 mm for an atom ve- 
locity of 815 m/s. The cavity spacing is adjusted by a 
piezoelectric modulator to be close to resonance with the 
1 So <-> 3 Pi transition by monitoring of a A = 791 .1 nm 
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probe beam through the cavity. The pump beam polar- 
ization direction is set perpendicular to the cavity axis 
to ensure maximum coupling. A magnetic field of w 150 
gauss oriented parallel to the pump beam polarization 
ensures that only m = <-> m = transitions occur. 

To measure the state of atoms which have passed 
through the pump beam, we monitor fluorescence of the 
atoms in the focused beam of a dye laser tuned to the 
1 So «-> Pi, A = 553.5 nm transition. The atoms were 
found to follow an adiabatic inversion process due to a 
Doppler frequency chirp from pump beam defocus. The 
effect is similar to that described in [9| . Incomplete inver- 
sion of the atoms by the pump beam is accounted for by 
defining an effective atom number N c g = N(p ee — p gg ), 
where p ee , p gg are the excited and ground state popu- 
lations of the atoms entering the cavity. For the data 
presented here p ee — p gg w 0.70. 

The atom-cavity interaction time ti n t = ^Ji:w m /vQ ~ 
0.10 fjs, with w m ps 41 pm the mode waist of the cavity 
and central atom velocity vo- 

In order to observe multiple thresholds, a well-defined 
atom-cavity interaction is essential, requiring a redesign 
of the earlier system 

Uniform atom-cavity interaction is achieved via control 
of atom-cavity coupling g and interaction time imt ■ The 
rapid spatial variation in coupling due to cavity stand- 
ing wave is eliminated by introducing a small tilt of the 
atomic beam relative to the normal to the cavity axis 
(cf. 53 )■ I n this configuration the microlaser exhibits 
two Doppler-shifted resonances approximately centered 
at cavity frequencies uj a ± kvoO where ui a is the atom 
resonance frequency. 

The peak traveling-wave atom-cavity coupling is given 
by g Q = (p/2h)^2Trhuj a /V w 190 kHz where V = 
ttLwq/4 is the cavity mode volume and p is the dipole 
matrix element. Note that this expression for go reflects 
a factor of 2 reduction in the peak coupling relative to 
the standing wave value. 

The atom beam is restricted to the center of the Gaus- 
sian mode by a 25/xm x 250/Ltm rectangular aperture lo- 
cated approximately 3 mm from the cavity axis. The 
aperture is translated and rotated in order to overlap 
the cavity axis as closely as possible. The resulting total 
atom-to-atom variation in peak coupling g is estimated 
to be fa 10%. 

To ensure uniform atom-cavity interaction times, a su- 
personic atomic beam with narrow velocity distribution 
is used. The beam oven is similar to that of Ref. 
and consists of a barium-filled resistively heated tanta- 
lum tube with a 250 micron diameter nozzle. The dis- 
tance between the oven and cavity is 45 cm. The longi- 
tudinal velocity distribution was measured via a Doppler 
fluorescence measurement to be maximum at vq ps 815 
m/s and have a width of Awfwhm ~ 100 m/s. A to- 
tal FWHM variation in gt- ln t of approximately 15% is 
achieved. Atomic density in the cavity is modulated by 
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FIG. 2: Gain (solid line) and loss (dashed line) per atom 
transit in the rate equation analysis (left and right sides of 
Eq. Q respectively) including effects of velocity distribution, 
nonuniform coupling, and detuning, with N c s = 1000. Closed 
circles: stable solutions. Open circles: unstable solutions. 



a translating knife edge approximately halfway between 
the oven and beam aperture. 

A heuristic rate equation for the microlaser may be 
obtained by equating the photon gain and loss per atom 
transit: 
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is the average emission probability for an atom and n is 
the number of photons in the cavity. The integrals over 
g and tint with respective weighting functions f(g) and 
h(ti n t) reflect the distributions in these values. (An anal- 
ogous stochastic averaging is justified in the appendix of 

The rate equation is expressed graphically in Fig. [21 
Damping in the oscillatory gain results from averaging 
over g and i m t. A solution n' of Eqn. ^ is stable if 
J7 (-Pcmit (n) ~ r c t int n/N cS )\ n , < 0. For sufficiently large 
atom number there exists more than one stable solution. 

We now describe two experiments to observe multiple 
thresholds, the transitions between solutions of Eqn. ^ 

In the first experiment, the microlaser cavity was 
locked on resonance (w C av = u> a + kvo9) and the mi- 
crolaser emission measured as a function of atom den- 
sity. The cavity was chopped between data collection and 
locking on resonance with the atoms. Locking was per- 
formed with a 791 nm probe beam tuned to uj a + kvoO via 
acousto-optic modulators and aligned with the TEMoo 
mode of the cavity. The cavity transmission was moni- 
tored by an avalanche photodiode (APD) and a feedback 
loop controlled the cavity PZT voltage to hold the cavity 
at resonance. Photons exiting the cavity were incident 
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FIG. 3: Photon number n vs. effective atom number N e g 
for cavity locking experiment. Circles (o): "unseeded" data; 
Crosses (x): "seeded" data. Solid line: Rate equation model, 
solution of Eqn.Qfor experimental parameters. 

on a cooled photomultiplier tube and pulses are recorded 
by a photon counter. 

We collected two sets of data: (i) in which the cavity 
photon number was allowed to develop from an initially 
empty cavity {n = 0), and (ii) in which the cavity was 
"seeded" with a large (n ~ 10 25 ) number of photons by 
a laser pulse provided by the cavity locking probe beam. 
In both cases the microlaser emission was observed for 
100 ms. 

Results from the cavity locking experiment are shown 
in Fig. [21 The solid line represents the solution to Eqn.Q] 
and is composed of two branches for the range of N e g 
shown. The first (lower) and second (upper) branches 
correspond to the first and second Rabi oscillations of 
Fig. [21 The lower half of the second branch (the portion 
of the line with negative slope in this figure) is unstable. 

Intracavity atom number was determined by measur- 
ing fluorescence at <-> 1 Pi, A = 553.5 nm with an 
imaging system and CCD camera. Uncertainties in cav- 
ity transmission, optical reflection coefficients, and CCD 
quantum efficiency make it difficult to measure the ab- 
solute density accurately. A fit to the experimental data 
was used to obtain a scale factor. A fitting parameter 
is close to the value expected from fluorescence measure- 
ments. The number of photons in the cavity is similarly 
estimated by considering the cavity parameters, losses in 
the optics, and detector quantum efficiency. In both cases 
the fitting parameters were within ±50% of expected val- 
ues. 

The unseeded results are in good agreement with the 



first solution branch. The seeded results agree with the 
first branch for densities up to the onset of bistability at 
N c g ~ 240; above this point the results agree with the 
second branch. 

A calculation for steady-state average photon number 
using single-atom micromaser theory of Filipowicz et al 
|l2| predicts a sharp transition between the first and sec- 
ond branches in at N® s w 415. No such transition was 
observed, apparently due to very long metastable life- 
times in the picture of the Fokker-Planck analysis of ^3] , 
due to the large photon number. 

In the second experiment, the cavity detuning A cav = 
w cav - u a is varied for constant atom density. Data for the 
detuning curves for a range of atom densities are shown in 
Fig. 01 The solid and dashed lines represent scans in the 
positive and negative detuning directions, respectively. 

Two-peaked structure is due to two traveling- wave res- 
onances of the tilted cavity. The Doppler splitting was 
measured to be about 27.1 MHz, corresponding to a cav- 
ity tilt angle 8 ss 13 mrad. 

As density is increased, the two resonances broaden 
and increase in amplitude. For N e g w 630 the second 
threshold appears as a "spike" near the peak of each res- 
onance, and at N e g sa 725 a jump to a third branch 
appears. 

The detuning lincshapes display strong asymmetry and 
hysteresis for atom numbers above the second thresh- 
old. For densities in which sudden jumps occur, each 
traveling-wave resonance of the cavity scanning line- 
shapes is highly asymmetric and shows hysteresis as 
a function of scan direction. This is most likely due 
to interactions with the two Doppler-shifted fields at 
uj a ± kvo9. For the general case in which neither field is 
dominant, the atoms experience a complex bichromatic 
interaction with these two fields. A simple model has 
been shown to qualitatively describe the lineshape asym- 
metry and hysteresis observed here. More detailed quan- 
titative descriptions are in progress. 

To examine the application of single-atom theory to 
our many-atom microlaser, a semiclassical theory |6j of 
the microlaser analogous to the Lamb theory of the con- 
ventional laser has been developed. In the case of 
monochromatic field and zero atom-cavity detuning, the 
semiclassical theory reduces to the rate equation model. 
This provides a more rigorous justification for Eqn. 

We have also performed quantum trajectory simula- 
tions of a microlaser with up to 5 intracavity atoms 
[I4j . The many-atom microlaser was shown to be in close 
agreement with the predictions from a single-atom the- 
ory |l2| , with a perturbation in the width of the photon 
number distribution due to cavity decay during the in- 
teraction time. 

The connection between the microlaser with the con- 
ventional laser may be illustrated by considering the ef- 
fect of increased variation in gt- mt . For simplicity suppose 
/(<?) = &{9~ 9o) an d interaction times are distributed ac- 
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FIG. 4: Cavity tuning lineshapes. Estimated photon num- 
ber n vs. atom-cavity detuning A cav . Solid lines, positive 
frequency scan (see text). Dashed lines, negative frequency 
scan. Estimated effective intracavity atoms: (a) N e g= 18 , 
(b) iV cff = 526, (c) N cff = 657, (d) N cS = 755. 



cording to atom lifetime 

■Pemit = TqT 1 J °° dt int exp 



The dimensionless gain is then 
t/r a ) sm 2 (y/rT- 

2„ 



l^int) = 

1/2 [((n + l) 5 2 r 2 )/(l + (n + l).g 2 r Q 2 )] . With Eqn.HJ this 
expression describes a conventional saturable gain laser, 
which is monostable and linear above threshold. 

In the microlaser, for any finite degree of broadening 
in gtint, we have P e mit{n — > oo) = 1/2. Therefore con- 



ventional laser behavior is recovered in the limit of large 
photon number. 

In summary, the microlaser is shown to exhibit multi- 
stability due to coherent atom-cavity interaction, which 
in conventional lasers is hidden by spontaneous emission 
and other incoherent processes. 

We have recently demonstrated that the microlaser 
field exhibits sub-Poisson photon statistics even for 
hundreds of intracavity atoms. Measurement of the spec- 
trum of microlaser emission 0] is also planned. 
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